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As confusion with lensing B-modes begins to limit experiments that search for primordial B-mode 
polarization, robust methods for delensing the CMB polarization sky are becoming increasingly 
important. We investigate in detail the possibility of delensing the CMB with the cosmic infrared 
background (CIB), emission from dusty star-forming galaxies that is an excellent tracer of the CMB 
lensing signal, in order to improve constraints on the tensor-to-scalar ratio r. We find that the maps 
of the CIB, such as current Planck satellite maps at 545 GHz, can be used to remove more than half 
of the lensing B-mode power. Calculating optimal combinations of different large-scale-structure 
tracers for delensing, we find that co-adding CIB data and external arcminute-resolution CMB 
lensing reconstruction can lead to significant additional improvements in delensing performance. 

We investigate whether measurement uncertainty in the CIB spectra will degrade the delensing 
performance if no model of the CIB spectra is assumed, and instead the CIB spectra are marginalized 
over, when constraining r. We find that such uncertainty does not significantly affect B-mode surveys 
smaller than a few thousand degrees. Even for larger surveys it causes only a moderate reduction 
in CIB delensing performance, especially if the surveys have high (arcminute) resolution, which 
allows self-calibration of the delensing procedure. Though further work on the impact of foreground 
residuals is required, our overall conclusions for delensing with current CIB data are optimistic: this 
delensing method can tighten constraints on r by a factor up to « 2.2, and by a factor up to « 4 
when combined with external « 3/iK-arcmin lensing reconstruction, without requiring the modeling 
of CIB properties. CIB delensing is thus a promising method for the upcoming generation of CMB 
polarization surveys. 


I. INTRODUCTION 

Many inflationary models make a key prediction: the 
production of a significant background of gravitational 
waves M. A detection of this stochastic gravitational 
wave background would therefore provide particularly 
strong evidence for the inflationary paradigm, along¬ 
side the observed super-horizon correlation, near scale- 
invariance and Gaussianity of the scalar density fluctua¬ 
tions (e.g., Q and references therein). Such a detection 
would also give significant insights into the physics of 
inflation. In particular, the ratio of gravitational wave 
or tensor power to scalar density power r gives a direct 
probe of the energy scale at which inflation occurred. 

The currently most promising method to measure this 
inflationary gravitational wave background is through 
its effects on the cosmic microwave background (CMB). 
However, tensor modes can be difficult to constrain ac¬ 
curately via the CMB temperature alone, because scalar 
modes also contribute to the temperature, which leads 
to cosmic variance and degeneracies. Unlike in temper¬ 
ature, the CMB sky contains two degrees of freedom in 
polarization (a magnitude of polarization and a direc¬ 
tion); because of this, one can construct a “null” linear 
combination from the two degrees of freedom which has 
no response to scalar perturbations at linear order, but 
still is sourced by inflationary gravitational waves. This 
combination, known as the CMB B-mode polarization. 
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therefore allows a constraint on primordial gravitational 
waves that is robust and free of leading-order scalar cos¬ 
mic variance 

Over the past year, the observational status of pri¬ 
mordial B-mode polarization has been a subject of ac¬ 
tive debate. While the BICEP2 collaboration reported 
a measurement of B-mode polarization on large angular 
scales 0, subsequent analyses based on Planck polar¬ 
ization measurements suggested that this signal could be 
sourced by foreground emission from polarized galactic 
dust A recently-released cross-correlation anal¬ 

ysis involving data from the BICEP2, Keck Array and 
Planck experiments [Tsjl has confirmed that a significant 
fraction of the signal seen in BICEP2 was of galactic ori¬ 
gin, and has thus replaced any evidence for inflationary 
gravitational waves with a new upper limit of r < 0.12 
at 95% confidence. For simplicity, and motivated by this 
new result, we work under the default assumption that 
r = 0 (though cosmic variance from tensor modes could 
straightforwardly be incorporated into our analysis if fu¬ 
ture measurements required this). While foreground po¬ 
larization appears to be a significant obstacle when ana¬ 
lyzing only one frequency channel, component separation 
with multifrequency data is a promising path towards re¬ 
moving much of the foreground contamination. 

Even in the absence of foregrounds, primordial gravi¬ 
tational waves are not the only source of B-mode polar¬ 
ization. A well-understood, recently observed [ill. [l5l - l2ll | 
source of B-mode polarization is the gravitational lensing 
of the CMB, the deflection of CMB photons by the dis¬ 
tribution of large-scale structure along the photons path. 
By remapping the CMB sky, gravitational lensing con- 
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verts some of the primordial E-mode polarization (the 
polarization component that can be sourced by scalars) 
into B-mode polarization . The lensing B-mode power 
spectrum is constant at low multipoles, so that the lens¬ 
ing B-modes induce variance and act as a source of noise 
similar in effect to instrumental white noise at a level of 
^ 4.9/xK-arcmin. The Keck Array / BICEP2 Collabo¬ 
ration recently reported the measurement of CMB maps 
that, for the first time, have noise levels below this value 
(3.4^K-arcmin, [l^). Lensing B-mode noise is thus be¬ 
ginning to limit even current B-mode surveys. Hence, 
for an improved constraint on r (or a better characteri¬ 
zation of the signal if it is non-zero), delensing methods 
are rapidly becoming essential [2J-[^. These methods 
involve subtracting off an estimate of the lensing B-mode 
map in order to reduce the lensing noise. 

To construct an estimate of the lensing B-mode and de¬ 
lens, one requires both a measurement of the E-mode po¬ 
larization and an estimate of the lensing which remapped 
the primordial CMB sky. Internal methods to delens in¬ 
volve reconstructing the lensing field (usiM statistical 
methods such as a quadratic estimator, [^, or a full 
maximum likelihood approach, [2^ [2^) with the same 
CMB dataset that is used to constrain the primordial B- 
modes. However, very low noise levels are required for the 
best internal delensing performance. Eurthermore, much 
experimental progress continues to be made with CMB 
experiments with large beams (^ 10'), which cannot re¬ 
construct CMB lensing maps that are sufficiently signal- 
dominated for efficient delensing. With such experiments 
reaching low noise levels, delensing which makes use of 
external datasets to estimate the lensing field may be 
promising. 

The cosmic infrared background (CIB) is diffuse far- 
infrared radiation originating in mostly unresolved dusty 
star-forming galaxies, where it is emitted by the UV- 
heated dust enshrouding young stars [1^. In a number 
of studies, the CIB was found to have a remarkabl y hig h 
correlation (si 80%) with CMB lensing (e.g., [13, l31|). 
confirming earlier predictions based on the similar high- 
redshift origin of both signals [Hj . This high correlation 
makes the CIB an ideal proxy for CMB lensing. In this 
paper, we will investigate the use of maps of the Cosmic 
Infrared Background to delens the B-mode sky. 

The use of external large-scale structure tracers for 
delensing was investigated in [13], focusing on tracers 
from surveys with well-known redshift origin. Eor such 
tracers, the authors concluded that external delensing 
was not a very promising technique, as futuristic surveys 
that mapped out all the structure to redshifts beyond 
z > 2 — 3 were required to achieve good delensing per¬ 
formance. Unlike the surveys considered, CIB maps are 
well-suited for delensing, because even currently avail¬ 
able CIB maps trace the underlying mass distribution 
over a broad range of high redshifts and thus are excel¬ 
lent lensing tracers [34l - l36l| . In a paper discussing how 
well delensing methods can improve constraints on the 
tensor tilt tit (which appeared as this paper was be¬ 


ing completed), [13, the authors included forecasts for 
a constant-correlation CIB channel. Our work differs be¬ 
cause it focuses on constraints on the tensor amplitude 
r rather than tensor tilt and assumes r is significantly 
less than 0.2. We also discuss additional aspects of CIB 
delensing, such as combinations of datasets, and, in par¬ 
ticular, the effects of uncertainties in the CIB cross- and 
auto-spectra. 

After briefly reviewing delensing in general, we discuss 
methods for CIB delensing and for calculating the 
resulting reduction in lensing noise in section 2. We then 
apply this formalism to models of the Planck CIB data in 
section 3 in order to forecast how much detection errors 
on r can be improved by CIB delensing. In section 3, we 
also discuss combinations of external delensing datasets: 
first, the optimal combination of CIB frequencies, and 
second, the combination of the CIB with external esti¬ 
mates of lensing from high-resolution CMB experiments. 
In section 4, we investigate the limitations of real data, 
focusing in particular on the problem of uncertainty 
in the CIB cross- and auto-spectra, and investigating 
to what extent CIB delensing can be self-calibrated 
without assuming models for the CIB spectra. We 
conclude in section 5. In this paper we assume a flat 
ACDM cosmology with flch^, Hq, lO^Ag, Ug, t} = 

{0.02226,0.1193,67.27,2.130,0.9653,0.063}, which is 
very close to the TT,TE,EE-blowP-|-lensing best-fit base 
ACDM model from [s^ - 


II. LENSING, DELENSING AND 
GORRELATION 

A. Preliminaries: Noise from Lensing B-modes 

As CMB photons travel through the universe, they 
are gravitationally lensed by large-scale structure. This 
lensing leads to a remapping of the CMB polarization 
anisotropies from their unlensed position n on the sky: 

l5(^) ~ Qunlensed(fl T d)j U (n) = Uunlensed(l^ T d) (1) 

with an analogous remapping in temperature. Here Q 
and U are the CMB polarization Stokes parameters, and 
d is the lensing deflection field. As the curl-like part 
of the deflection is many orders of magnitude smaller 
than the gradient-like part [13 , the lensing can be equiv¬ 
alently described using the convergence k = — • d, 

which can be shown to equal a projection of the matter 
density out to high redshift [13 ■ The convergence is also 
simply related to the often-used lensing potential (j) by 

Writing the polarization not in terms of the experimen¬ 
tally measured Q and U parameters but in terms of even 
parity E and odd parity B-modes, we find that lensing 
induces a conversion of a pure E mode pattern into B- 
modes, with the lensing B-mode given in harmonic space 
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by (to leading order in k): 

= I wvHi - 1 ') ( 2 ) 

for 

21'. n — I'i 

W^(M') = ^p3^sin(2(^i.,0, (3) 

where E is the unlensed E-mode, /, 1 are multipoles writ¬ 
ten in the flat sky approximation, and (pi y is the angle 
between 1 and 1'. The flat-sky approximation used here 
gives sufficient accuracy for our simple forecast calcu¬ 
lations and Trange, though we also verify our numeri¬ 
cal results using the corresponding full-sky formulae [i^ ■ 
Keeping only the lowest order term in n has also been 
found to be an excellent ap proximation on the scales rel¬ 
evant for delensing [H, |4l| . 

From this expression, we can calculate the power spec¬ 
trum 


(^lens(i)^le„s* ^ ( 27 r) 2 ^^(l - 


to give 


J ^ ^2 ^ ( 5 ) 

This lensing B-mode power is effectively constant on large 
scales Z < 100 [s^. 

The full (beam-deconvolved) B-mode power spectrum 
is given by 

^ ^BB.r ^ ^SB.lens j^BB 

where is the contribution of tensor modes to the 

B-mode power spectrum (assumed zero by default) and 
is the instrumental noise power 

= (Ap/Tcmb)^ e'"®PWHM/( 81 n 2 ) 

with Ap the noise level, here written in /iK-radians 
(though typically specified in /xK-arcmin), and 0fwhm 
the full width at half maximum of the beam in radians. 
Assuming the realization of the lensing B-modes on the 
sky is unknown, the presence of lensing B-mode power 
will contribute to the error on a measurement of the B- 
mode power spectrum: 



where /sky gives the sky fraction observed by the sur¬ 
vey. While non-Gaussian covariance between different B- 
multipoles is present because fluctuations in lensing and 
E-modes simultaneously affect different B-mode scales 
this covariance is subdominant on the large 
scales that constrain r. Despite this, we include this co- 
variance in our analysis in later sections of this paper. 


The error on the B-mode power determines the error on 
a constraint of the tensor-scalar ratio r (assuming the 
“detection” null hypothesis of no tensor modes in the 
error calculation): 
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In the second line we have assumed that both the lens¬ 
ing B-mode and noise power spectra are constant, i.e., 
white-noise-like, at I < 100 [s^, where nearly all the con¬ 
straining power on r arises (note that the derivative with 
respect to r is simple since r merely normalizes a fixed 
tensor spectrum). The lensing-induced error affects any 
measurement of the B-mode power spectrum and there¬ 
fore any constraint on r, and persists even as instrumen¬ 
tal noise is reduced to zero. Delensing algorithms are 
therefore needed to remove the lensing B-mode and thus 
reduce the errors when the instrumental noise approaches 
or falls below the level of (jBB.iens^ which corresponds to 
Ri 4.9/iK-arcmin in polarization. 


r(r) = 


dC 


BB,7 
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B. Delensing with Large-Scale Structure: Methods 
and Expected Performance 

We now consider delensing with a cosmic infrared back¬ 
ground map /, which traces the underlying matter distri¬ 
bution and which we will attempt to use as a proxy for k. 
(Our analysis is general and thus I could be any tracer of 
large-scale structure.) We construct an estimate of the 
lensing B-mode using this field, using a weighting /(1,1') 
(which we will determine): 

!')/(!, - 1 ') ( 10 ) 

where E^ is the noisy, observed E-mode; to lowest order, 
the fact that it is lensed can be neglected. Though we 
distinguish E and E^, we will always take I to contain 
all the components of the observed GIB map (e.g., shot 
noise, foregrounds, instrumental noise if relevant). 

We then subtract this estimate off from the measured 
B-mode in order to delens. The residual lensing B-mode 
is given by 

- f 

^^"(l) = B'™"(1) - b1“"(1) = J X 

{E{1')k{1 - 1') - /(I, l')E^(l')/(l - 1')) (11) 

We now wish to determine / such that the residual 
lensing B-mode power, and hence the lensing B-mode 
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noise, is minimized. Evaluating the lensing B-mode 
power spectrum after the delensing term has been sub¬ 
tracted gives: 

^BB,Tes _ / ^ ^ 

-(/(i,i') + r(i,i'))Q^^qiiFi (12) 

+r (1, i')/(i, 

Taking the derivative with respect to /, we obtain the 
form of / which minimizes the residual lensing B-mode: 








■Np 


EE 




(13) 


We note that for signal-dominated E-modes, this weight¬ 
ing corresponds simply to a normalization of the I field 


by Cl /Cl . Additionally, noisy E-modes are down¬ 
weighted. We also note that the same weight / could be 
obtained by treating the part of the GIB that is uncorre¬ 
lated with lensing as effective noise in the standard cal¬ 
culation of internal delensing methods as in 331 and [3^ . 
However, our explicit approach has advantages when we 
consider uncertainty in the GIB spectra in later sections 
of this paper. 

Using this (perturbatively) optimal weight /, we ob¬ 
tain the following residual lensing B-mode power after 
delensing: 


C, 


BB,res 

I 


(14) 


1 - 




Cj/E + ^EE ) P\l-V\^ 


where p is the correlation coefficient of GMB lensing and 
the GIB: 


Pi 


cf 

^crc/^ 


(15) 


The full delensed B-mode power spectrum then becomes 
(from Eq. 6): 

Cf ®’dei = Cf -b Cf -b NP^. (16) 

When the E-mode is signal dominated for all relevant V 
used in the sum in Eq. 14, the delensing can be simply 
described by the following reduction of the lensing signal 
in the map: 


cr ^ (1 - ppcr- ( 17 ) 

If the correlation coefficient is constant with I , the lensing 
B-mode power and the resulting lensing noise are then 
similarly reduced by a factor (1 — p^). 


III. DELENSING WITH THE GIB 

A. Modeling the CIB-Lensing Correlation for 
Forecasting Delensing Performance 

To estimate the delensing performance and the reduc¬ 
tion in noise for B-mode searches, we must know the 
residual B-mode power; to determine the residual B- 
mode power with Eq. 14, we must in turn have an es¬ 
timate for the correlation coefficient of the GIB with 
GMB lensing, pi. To calculate pi in this work, we will 
use the GIB halo model of to estimate GIB lensing- 
cross- and auto-power spectra; this model is a good fit 
to Planck measurements of these spectra [13,143. As the 
GIB power spectra in the default model do not include 
the shot noise contributions that remain in the measured 
spectra, we add shot noise terms for both bright GIB and 
radio sources; we obtain the relevant amplitudes (which 
match the measured spectra) from the models used in 
[ 43 . However, in this section we do not include instru¬ 
mental noise (or any additional noise due to component 
separation). On the scales relevant for delensing, this 
is an excellent approximation for the higher frequency 
Planck channels (545, 857 GHz), but not necessarily the 
lower frequency channels (353 GHz, 217 GHz), which are 
only included for comparison. We will discuss the impact 
of uncertainties of GIB spectra in section 4. 

The correlation coefficients for the different frequencies 
are shown in Fig. 1. As expected from prior work [30| . 
we find a high degree of correlation at low multipoles, 
peaking above 0.8 for all the frequencies. This again 
illustrates that the GIB is an excellent tracer of lensing. 
The correlation coefficients do not a ppe ar constant with 
scale, unlike shown in Figure 13 of |30|; in part, this is 
because we include the shot noise contribution to the 
GIB power spectra, though even without the inclusion 
of shot noise the correlation in the model we use falls 
somewhat with 1. We will assume that only GIB scales 
I > 60 can be used to delens; however, the results are 
not very sensitive to the exact value of the \ow-l cutoff, 
as most of the lensing B-modes arise from higher-Z lens 
scales (see Fig. 1). 

We now consider to what extent foregrounds, in par¬ 
ticular galactic dust emission, will reduce the correlation 
from the models described above and degrade the de¬ 
lensing performance. To determine this, we consider the 
effect of foregrounds on the GIB maps’ cross-correlations 
with GMB lensing as well as on their power spectra. 
Galactic foregrounds are not expected to be correlated 
with the true lensing maps (so that the cross-power mea¬ 
sured in [lO] should be accurate everywhere), but they 
will contribute to the GIB auto-power spectra on large 
scales. However, on « 2240 deg^ of sky, the Planck 
Gollaboration have demonstrated that galactic dust 
can be removed by subtracting dust maps constructed 
from observations of galactic HI emission (21cm emis¬ 
sion from neutral hydrogen, [13) combined with Planck 
HFI and IRAS observations; after subtraction, the dust 
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FIG. 1. Correlation coefficients for the different Planck CIB 
frequencies (545, 857, 217, 353 GHz) with GMB lensing, us¬ 
ing the best-fit model from [^. We also show the correlation 
coefficient for the co-added linear combination of these CIB 
frequencies. The contribution of shot noise to the correlation 
coefficient is included (at a level which matches the measured 
CIB spectra). We neglect instrumental noise (or any addi¬ 
tional noise due to component separation), which is an ex¬ 
cellent approximation for the high frequency Planck channels 
(545, 857 GHz), but not necessarily the lower frequency chan¬ 
nels; these are thus only included for comparison. Finally, 
plotted beneath the correlation curves, the black solid line 
shows (Gf*^ X 9G®®/9Gf'‘);/<ioo (in arbitrary units), i.e. how 
much each scale I in the lensing contributes to the mean B- 
mode power on scales I' < 100, and thus which scales are most 
important to delens. 


contributions to the auto-spectra are small on angular 
scales above I = 150. While the Planck CIB analy¬ 
ses that demonstrate galactic dust subtraction rely on 
some of the cleanest areas of sky, most upcoming B-mode 
searches will either focus on these same low-dust regions 
at high galactic latitude where foreground contamination 
is smallest, or will have enhanced multifrequency cover¬ 
age that will allow improved dust modeling and subtrac¬ 
tion. Even if large scales I < 150 are assumed completely 
unusable, the reduction in delensing performance when 
changing from = 60 to ^min = 150 is only si 12%, 
because smaller scales I = 200 — 500, which are less af¬ 
fected by dust, are much more important for delensing 
(as can be seen in Fig. 1). Though for the largest surveys 
our forecasts may be somewhat optimistic, our correla¬ 
tion model should thus be sufficiently accurate to make 
realistic forecasts of delensing performance for most sur¬ 
veys. In our analysis we have assumed that the E and 
B GMB polarization maps have been cleaned with mul¬ 
tifrequency data so the levels of residual galactic dust 
and CIB contamination in the CMB maps are negligi¬ 
ble, though we revisit this in section 4, along with other 


Frequency 

Error Improvement Factor a 

857 GHz 

2.2 

545 GHz 

2.2 

353 GHz 

2.1 

217 GHz 

2.0 

545 -t 857 GHz 

2.3 

all frequencies coadded 

2.3 


TABLE I. Ideal improvement factors on errors on a detec¬ 
tion of r achievable with CIB delensing using four Planck 
frequency channels, assuming zero instrumental noise in both 
the CMB and the CIB maps (a good approximation for 545 
and 857 GHz, but not for the lower frequencies). 

aspects of potential foreground contamination. 

B. Delensing Performance: Individual Frequencies 

With this model for the correlation coefficients, we can 
characterize the delensing performance for CIB maps at 
different frequencies. We calculate the delensed resid¬ 
ual B-mode power for the different Planck CIB channels 
from Eq. 14, at first setting the instrumental noise in the 
CMB maps to zero. The residual lensing B-mode spec¬ 
tra, after delensing in this idealized case, are shown in 
Fig. 2; we find that for all Planck frequencies CIB de¬ 
lensing can remove a substantial part of the lensing-B 
mode power. We quantify this by measuring the change 
in the mean delensing, averaged over large 

scales I < 100, where most of the constraining power 
for primordial B-mode polarization lies. We obtain the 
following results: after delensing noise-free CMB maps 
using the 857 GHz, 545 GHz, 353 GHz, and 217 GHz 
GIB frequencies, the mean low-ell reduced 

by a factor 0.46, 0.45, 0.47, 0.49 respectively (though the 
217 and 353 GHz delensing performance will be degraded 
when noise is included.) With our knowledge of the resid¬ 
ual B-mode power, we can now define an improvement 
factor in the measurement error of r, a{r), as in [3^ : 

a = Cro(r)/CTdelensed(?') (18) 

where a{r) depends on the B-mode power as in Eq. 9, uo 
indicates the error before delensing, and (Tdeiensed is the 
error after delensing (i.e., determined by rather 

than When the instrumental noise in Eq. 9 

is neglected, this reduces to the simple expression a = 
^^BS,iens)^^i^o/((7^^^>"'®®)j<ioo because both the lensed 
and residual B-mode power are flat at I < 100 (see Fig. 2); 
we show the results for this expression (which is just the 
inverse of the reduction factors quoted earlier) in Table 
1 . 

Going beyond the ideal noise-free limit, we now ex¬ 
amine the improvement in measurement errors on r as 
a function of instrumental noise level of a CMB polar- 
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ization telescope. With non-zero experimental noise, the 
improvement factor becomes: 




BB 


[Ap]) 


/<100 




(19) 


where the residual B-mode power is calculated from 
Eq. 14 as before, now including the effects of instrumen¬ 
tal noise in the calculation. 

Since the noise power spectrum after beam deconvolu¬ 
tion enters into this calculation, the impact of noise de¬ 
pends on the instrumental beam size (as seen in Eq. 7). 
We will thus consider two reference experiments in our 
analysis. The first reference experiment, which we will 
designate the highRes experiment, has high resolution - a 
1.4 arcmin FWHM Gaussian beam - and is thus capable 
of efficient internal delensing when noise levels are low. 
This experiment should be taken as a proxy for ground- 
based telescopes such as ACTPoL POLARBEAR-1/11, 
SPTpol, and their successors Eiill , as well as high- 
resolution satellites such as COrE or PRISM Is^ . 
The second reference experiment, which we term the 
lowRes experiment, is a lower angular resolution CMB 
telescope targeted towards measuring large-scale CMB 
polarization. We take it to have a 30 arcmin FWHM 
Gaussian beam; this resolution is too low for efficient in¬ 
ternal delensing. The lowRes experiment is intended to 
represent CMB telescopes such as Keck Array, BICEPS 
or LiteBIRD For the lowRes experiment, we as¬ 

sume that the measured E-mode can be co-added with 
Planck maps of the E-mode polarization on smaller scales 
(for which we assume bOjuK-arcmin in polarization and a 
7-arcmin beam, as in |57l|l. 

For both experiments, we show the improvement fac¬ 
tor a for delensing with the 545 GHz CIB frequency in 
Fig. 3, as a function of polarization noise level. Other 
CIB frequencies give nearly the same results, as their 
correlation coefficients with leasing are very similar. For 
the highRes experiment, we see that a becomes greater 
than 1.5 as noise levels fall below 4/rK-arcmin and reaches 
Ri 2.2 as the experimental noise approaches zero. As a 
CIB-delensed survey gives the same constraints as an 
times larger survey that has not been delensed (neglect¬ 
ing windowing effects), these improvement factors can 
correspond to signihcant enhancements in experimental 
power. For the lowRes experiment, the delensing perfor¬ 
mance, though similar, is slightly reduced for intermedi¬ 
ate noise levels. This is due to differences in the E-mode 
map: while for the highRes experiment, the E-mode map 
is signal-dominated on all relevant scales (i.e., the results 
do not change visibly if we set the noise in the E-mode 
map to zero), the lowRes E-mode is not signal-dominated 
everywhere, which results in a small reduction in delens¬ 
ing performance. 

For the highRes experiment, we also show the improve¬ 
ment factor a for internal delensing, i.e. for delensing 
with an estimate of the leasing convergence reconstructed 
from CMB data. We calculate the expected noise level 
of CMB leasing reconstruction for the minimum variance 


combination of optimal quadratic estimators as in 
(assuming scales of I = 100 — 3000 and I — 100 — 4000 are 
used in temperature and polarization reconstruction re¬ 
spectively); we hence obtain a correlation coefficient for 
the reconstructed maps, which allows us to calculate the 
improvement factor a as before and compare with CIB 
delensing. 

In Fig. 3, we see that for noise levels above Ri 5 — 6/rK- 
arcmin, CIB delensing performs slightly better (though 
a is not large); for much lower noise levels, with a small 
beam such that lensing reconstruction is possible, in¬ 
ternal leasing reconstruction performs better than the 
CIB at delensing the highRes experiment. The use of an 
iterative procedure for internal reconstruction will not 
greatly modify the reconstruction curves at the noise 
levels shown in Fig. 3 [33|, though it will cause large 
increases to a when the noise is below Ri 2/iK-arcmin. 
It is encouraging that for moderate noise levels > 2fiK- 
arcmin, internal delensing and CIB delensing (with cur¬ 
rent CIB data) are fairly similar in performance. Fur¬ 
thermore, for future lowRes-\ike experiments without the 
necessary angular resolution to perform lensing recon¬ 
struction, CIB delensing will allow an improvement at 
low noise levels. 



FIG. 2. B-mode power spectra after delensing with CIB maps 
at different frequencies (and with a linear combination of 
maps at all frequencies). The upper solid line shows the orig¬ 
inal lensing B-mode power spectrum. We here assume that 
the CMB maps have negligible instrumental noise on all rel¬ 
evant scales. As in Fig. 1, we also neglect instrumental and 
component separation noise in the CIB maps (which is a good 
approximation for the Planck 545 and 857 GHz maps). 


C. Combining Different CIB Frequencies 

We now investigate how well we can delens using com¬ 
binations of large-scale-structure datasets, turning first 
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Ap [/^K—arcmin] 

FIG. 3. Improvement factors in errors on a detection of r 
when delensing with the 545 GHz GIB channel, as a function 
of the polarization noise level of the experiment. The blue 
solid line shows GIB delensing performance for an experiment 
with a 1.4 arcmin beam (our highRes reference experiment). 
The green dashed line shows results for an experiment with 
a 30 arcmin beam (the lowRes experiment). For this exper¬ 
iment, we assume the E-mode map used for delensing is ob¬ 
tained in combination with E-mode maps from Planck. It can 
be seen that for low noise levels, GIB delensing can improve 
constraints on r by a factor of up to « 2.2. For comparison, 
we also plot the improvement when delensing the highRes ex¬ 
periment using internally reconstructed (minimum variance 
quadratic estimator) lensing maps (this requires high resolu¬ 
tion). a is shown both for lensing reconstruction alone (dotted 
black line), and for lens reconstruction coadded with the 545 
GHz GIB (dotted red line). 


to combinations of different GIB channels. Writing our 
combined delensing tracer / as a linear combination 
/ = 'Y^Cili, where the fields li are different GIB fre¬ 
quencies (or other large-scale structure tracers), we must 
determine the coefficients Ci that allow us to delens op¬ 
timally. We note that in Eq. 14, these coefficients only 
enter into the correlation coefficient of I with lensing, pi; 
hence, if we seek to maximize the delensing performance, 
we must maximize pi with respect to the Ci coefficients. 
In the appendix, we solve for the coefficients Ci that max¬ 
imize the correlation. We find that for the optimal linear 
combination I — ^Cilt, 




jPjK.\ 


^K.H 


C, 


hh 


( 20 ) 


where p*” is the cross-correlation coefficient of large-scale 
structure tracer li with the true lensing convergence, 
is the power spectrum of li, and p^^ are the cross 
correlation coefficients of tracers li and Ij. 

For the linear combination I, the correlation coefficient 


with GMB lensing is 

P — ^ ^ PiniP )ijPjK.- (21) 

ij 


at each 1. 

For the linear combination of different GIB frequen¬ 
cies, we have all the information needed to calculate this 
correlation coefficient from the model for the GIB we use 
in this paper. We consider the combination of 545, 857, 
217 and 353 GHz GIB maps, as well as the combination 
of only 545 and 857 GHz maps; the cross-correlation co¬ 
efficient for the combination of four frequencies is shown 
in Fig. 1. For both linear combinations, we note that 
the improvement in correlation is quite small compared 
to individual GIB frequencies. This is expected: given 
how correlated the different GIB frequency channels we 
consider are (typically ~ 95%, 113), it is clear they must 
contain much of the same information. 

Nevertheless, we can calculate the delensing improve¬ 
ment factor a for these linear combinations. As for indi¬ 
vidual frequencies, we obtain the improvement factors a 
by calculating the delensed residual B-mode power from 
Eq. 14. As can be seen in Table 1 and as expected from 
the correlation coefficient, the value of a for a noise-free 
experiment is 2.3 for both 545 GHz -k 857 GHz and for 
all GIB frequencies combined (neglecting any GIB noise) 
- not very different from the results of delensing with the 
545 GHz channel alone. 

Larger improvements may be possible if a wider range 
of less correlated frequencies is used, but we defer an 
investigation to future work. 


D. GIB Maps Combined with Lensing 
Reconstruction 


We now consider combining maps of the GIB with 
maps of lensing that have been reconstructed from GMB 
data. Such a combination could provide a better proxy 
for the true GMB lensing field and thus increase the 
delensing performance. To derive the delensing perfor¬ 
mance of such a combination, let us assume that our 
lensing reconstruction, for a certain instrumental noise 
level and beam, results in a lensing map with reconstruc¬ 
tion noise , so that its power spectrum is 

(^KrecKrec _ (JKK. J^KK. ^ 22 ) 

We combine this reconstructed lensing map with a GIB 
map at one frequency (e.g., at 545 GHz), which we take 
to have a correlation coefficient p\ with the true lensing 
field. How efficient is this combination at delensing? To 
calculate this, we treat our lensing reconstruction as just 
another large-scale structure tracer A and repeat our cal¬ 
culation of how well a combination of tracers can delens. 
The correlation matrices of Eq. 21 can be simply calcu- 









lated to give the total correlation coefficient (at each 1): 


(1-(P?)^)§- 


Pi 


= 1 - 


(1-(P?)") + 


Nr 


(23) 


Given our knowledge of the CIB-lensing correlation co¬ 
efficient and our knowledge of the lensing reconstruction 
noise, we can calculate the improvement factor a. 

We will first consider the case where the CMB lens¬ 
ing is estimated internally from the same survey that 
is being delensed. This requires an experiment with a 
small beam, such as our highRes reference experiment. 
In Fig. 3, we include a curve showing the improvement 
factor a for delensing highRes with the combination of 
545 GHz GIB maps and quadratic estimator lensing re¬ 
construction maps; moderate increases in a can be seen 
when delensing with this combination. 

We now consider the delensing synergy of two differ¬ 
ent overlapping surveys - an experiment similar to the 
lowRes experiment (with low noise and low resolution 
that constrains large-scale B-modes and r, but cannot 
reconstruct lensing well) and an experiment similar to 
the highRes experiment (with high resolution that recon¬ 
structs K as a delensing tracer). This is a useful scenario 
because the high resolution experiment can be used to 
delens the low-resolution experiment if the surveys over¬ 
lap [1^. The polarization noise level of the high reso¬ 
lution experiment /\^*shRes course generally differ 

from that of the low resolution experiment We 

show the delensing performance for such a scenario as 
a function of in Fig. 4 (for highRes noise val¬ 

ues of Ap®^^®'* = 2.8/iK-arcmin and Ap®*'^®® = 7.1/rK- 
arcmin). The lensing reconstruction from highRes is ei¬ 
ther used alone or combined with the 545 GHz GIB maps 
to delens lowRes. 

We find that the combination of GIB and reconstructed 
lensing from high-resolution GMB maps, even with inter¬ 
mediate noise levels of Ap®*'^®® si 3/iK-arcmin, is quite 
effective as a delensing tracer, and is able to improve 
the detection errors by close to a factor 4 for a lensing- 
limited survey. The addition of the GIB results in a sig¬ 
nificant improvement in delensing performance compared 
to reconstruction with the highRes experiment alone; 
noisy Ap®*'^®® = 7.1/iK-arcmin reconstruction coadded 
with GIB data delenses as well as Ap®*'^®® = 2.8/iK- 
arcmin reconstruction without the GIB. This increased 
performance is due to the fact that lensing reconstruc¬ 
tion at intermediate noise levels does not image small 
lenses well; adding the GIB provides information on these 
small lenses and increases the small-scale correlation of 
the coadded delensing tracer with the true lensing con¬ 
vergence. The combination of a low noise B-mode search 
(similar to the lowRes experiment) that cannot recon¬ 
struct lensing well with a higher resolution experiment 
(such as the highRes experiment), even at a higher noise 
level, is thus quite promising for delensing, especially 
when this external lensing reconstruction is combined 



^bwRes J^x_a,rcmin] 


FIG. 4. Delensing improvement factor a for delensing a 
low-noise/low-resolution CMB polarization experiment (e.g.. 
Keck Array/BICEP3 - for us, the lowRes reference ex¬ 
periment) with lensing maps reconstructed from an over¬ 
lapping higher resolution CMB experiment (e.g., ACT- 
Pol/SPTpol/POLARBEARI/II - for us, the highRes refer¬ 
ence experiment). The lensing reconstruction maps are used 
either alone or in combination with 545 GHz CIB maps. We 
assume all constraints on r are from the lowRes experiment, 
which has a noise level given by Ap'"^®°. The external CMB 
lensing reconstruction used for delensing is assumed to be 
reconstructed with a quadratic estimator from the highRes 
experiment, which has a noise level Shown are 

values of Ap®^^®® = 2.8/tK-arcmin (upper dashed line) or 
Apgh^e® = 7.1/iK-arcmin (lower dashed line). The solid lines 
correspond to delensing with optimal combinations of lens¬ 
ing reconstruction and the 545 GHz Planck CIB map (with 
the same ordering by highRes noise). It can be seen that 
adding CIB data to a lensing reconstruction map greatly in¬ 
creases the delensing performance, and that delensing a low- 
resolution, low-noise experiment with a combination of CIB 
maps and overlapping intermediate-noise external lensing re¬ 
construction is promising. We assume that the E-mode map 
used for delensing is signal dominated, which is likely to be 
the case if we can combine high- and low-resolution data. 


with GIB data. (We note that our calculations are only 
relevant when the constraint on r from the lowRes ex¬ 
periment is more powerful than the constraint from the 
highRes experiment alone, e.g., because the latter has 
higher noise or has not been optimized to measure large- 
scale polarization). 

Additional improvements might be obtained by com¬ 
bining GIB data with other tracers of large-scale struc¬ 
ture, such as galaxy surveys (or other GIB maps over a 
wider frequency range), in order to improve the overlap 
of the combination’s redshift origin with the lensing ker¬ 
nel. However, we leave such an investigation to future 
work. 








9 


IV. POTENTIAL LIMITATIONS OF 
OBSERVATIONAL GIB DATA AS A DELENSING 
TRACER 


A. Uncertainty in the GIB Cross- and Auto-Power 
Spectra, and Prospects for Self-Calibration 


In our calculations of section 3, we assumed a model 
(albeit one based on measurements) for the lensing cross- 
and auto-spectra and of the CIB. However, 
though good measurements are available, the true cross- 
and auto-spectra of the CIB fields are of course not ex¬ 
actly known. This is potentially important for estimating 
delensing performance, not just because forecasts may 
be slightly inaccurate, but because the true CIB lensing- 
cross- and auto-spectra determine the residual B-mode 
power spectrum (the power spectrum of the lensing B- 
mode that remains after delensing) as in Eq. 12. Be¬ 
cause the CIB spectra affect the residual B-mode power 
spectrum, insufficiently constrained CIB spectra can de¬ 
grade constraints on r. For example, if the CIB cross¬ 
correlation is in fact smaller than in the fiducial model, 
the reduced correlation leads to an “excess” of residual 
lensing B-mode power compared to the fiducial model, 
which could be confused with the presence of primordial 
tensor B-modes. We illustrate this by constructing 10 dif¬ 
ferent CIB cross-correlation curves, all of which are con¬ 
sistent with current Planck measurements, and calculat¬ 
ing the residual B-mode power spectrum from Eq. 12 for 
each such model. (To construct each curve, we perturb 
every CIB-lensing-cross model bandpower by a random 
amount drawn from a Gaussian of width given roughly by 
the size of the measurement errors from (30j ; see later for 
more discussion of our error assumptions). In Fig. 5 we 
show the difference between the residual B-mode spec¬ 
trum for each of these 10 alternative “models” and the 
residual B-mode spectrum for the fiducial model. We also 
show a primordial tensor B-mode power spectrum corre¬ 
sponding to r = 0.001 for comparison. Though the dif¬ 
ferences in the residual B-mode power are small (of order 
5% of the total power), it can be seen that they are com¬ 
parable in size to primordial B-modes when r ~ 0.001; 
uncertainty in the CIB spectra could thus limit the power 
of CIB delensing for B-mode surveys that are capable of 
such high precision constraints on r (e.g., are sufficiently 
large). 

We now seek to quantify the effects of uncertainty in 
the CIB spectra on the delensing procedure more rigor¬ 
ously. We begin by writing a theoretical model for the 
full delensed B-mode power spectrum, as: 


Cf a,,h) = (a,), C"(6,)] 


(24) 

where we have denoted the residual lensing B-mode 
power as a function of the CIB spectra, and we have 
assumed a measurement without noise bias. We have 
also assumed that, in turn, the parameters ai,bi (which 
we will later define as effectively parametrizing the am¬ 



FIG. 5. Illustration of the possible impact of uncertainty 
in the true CIB-lensing cross-correlation on the residual (de¬ 
lensed) B-mode power spectrum. We construct 10 different 
possible “models” of the true cross-power by perturbing band- 
powers of our fiducial model by a random amount drawn from 
a Gaussian, Planck-like error distribution. Each model cross¬ 
power is consistent with current measurements but results in 
a different residual B-mode power spectrum; the excess resid¬ 
ual B-mode power compared with our fiducial model is shown 
in blue. For comparison, the power spectrum of primordial 
B-modes with an amplitude corresponding to r = 0.001 is 
shown. It can be seen that variation in the uncertain CIB 
spectra can be degenerate with the presence of primordial 
B-modes, especially if only the largest scales of the B-mode 
polarization are measured. 


plitude of CIB bandpowers) completely describe the true 
CIB cross- and auto-spectra. As illustrated previously, 
uncertain CIB spectral parameters Ui, bi could poten¬ 
tially be degenerate with r and thus reduce the perfor¬ 
mance of CIB delensing. To quantify this limitation, we 
will again forecast constraints on r, but now marginal¬ 
ize them over the uncertain true CIB spectra (i.e., the 
parameters ai,bi). 

First, we discuss how the CIB spectra can be 
parametrized for delensing analyses. While considerable 
progress has been made in modeling these spectra, for 
cosmological constraints on r we believe it is more robust 
if cross- and auto-spectra can be completely constrained 
from measurements without assuming a certain model for 
the CIB properties and spectra. This is especially impor¬ 
tant as foreground residuals may also contribute to the 
spectra and may not be simple to model. We thus pro¬ 
pose to constrain r and the CIB spectral parameters ai, bi 
simultaneously from measurements of as well as 

measurements of the CIB spectra and and then 
marginalize constraints on r over ai,bi. should ide¬ 
ally be measured with the exact CIB map used for delens¬ 
ing, including any foreground (e.g., galactic dust) residu- 
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als present in the map. As long as fiducial model spectra 
are close to the data, they should be used for the filtering 
function / of Eq. 12. This preserves linearity and is only 
slightly suboptimal, because any correction to the de- 
lensing performance is second order in an error in /. For 
such an empirical approach, we must have a very general 
parametrization of the CIB spectra. As indicated pre¬ 
viously, we will here assume that the CIB lensing-cross- 
and auto powers can be completely described by rescal¬ 
ing different /-ranges of the fiducial CIB spectra 
similar to rescaling effective bandpowers: 

Cf{a,) = Cf{l + a,Um, 

Cl\h) = Cl^ {l + hUm (25) 

where a sum over i is implicit, and Ui are top-hat func¬ 
tions that are unity within an /-range corresponding to 
the ith bin (we describe our binning in later sections), 
and zero otherwise, i.e.: 


[/,(/) 


1 

0 


If ^i,min ^ / ^ /i,max7 

otherwise. 


(26) 


(though there are other possible, equally complete, 
parametrizations). In writing this model, we assume that 
the true CIB spectra in our data are sufficiently smooth 
that they can be described without significant losses by 
this discrete parametrization when integrated over the 
kernel of Eq. 5 (and Fig. I), though this assumption ap¬ 
pears reasonable unless the binning is very coarse. 

Assuming that measurements of both CIB spectra and 
the B-mode power spectrum are available, and that we 
are fitting both r and the CIB spectral parameters simul¬ 
taneously, we now investigate to what extent uncertain 
Ci^ can degrade the constraints on r. We will neglect un¬ 
certainty in the CIB auto-power, because we assume that 
the raw auto-spectrum of the CIB map (including any 
residual foregrounds) can be measured with Planck data 
much more accurately than the leasing cross-spectrum. 
This is currently the case for the relevant high frequency 
CIB maps at 545 and 857 GHz (despite the smaller area 
on which the auto-spectrum is measured), and is typi¬ 
cally a good approximation because the leasing recon¬ 
struction maps are expected to have significantly higher 
noise than the CIB maps. (Furthermore, the residual 
B-mode is twice as sensitive to changes in the cross¬ 
power than to the auto-power, as can be inferred from 
Eq. 12). To investigate how much constraints on r are 
degraded by uncertain , we perform a Fisher forecast 
for the parameters 9p = (r, a^), obtaining a constraint on 
r marginalized over the CIB-lensing cross-power param¬ 
eters Ui- We construct the Fisher matrix as 


,BB 


E E 


da 


SB,del 


+ E 


mm ^ min 

dCf dCf 
'Wr~d9^ 

(ACfy 


de. 


[C' 


OV 




-1 dCl 


BB,del 




dOa 


where the cross-power has been binned into bandpow¬ 
ers j which are identical to the effective bins of the ai 
parametrization. To construct this Fisher matrix, we 
must evaluate derivatives of and with re¬ 

spect to all parameters Op, and estimate the measure¬ 
ment errors (covariance matrices) of the measurements 
of Cf and CfC 

We begin by calculating how changing the assumed 
CIB cross-spectrum via Ui affects the delensed residual 
B-mode power. Recalling our expression for the residual 
B-mode power (Eq. 12), we can insert our parametriza¬ 
tion of the CIB spectra (Eq. 25) in terms of ai,bi into 
this expression (keeping the weight / fixed to the fidu¬ 
cial model) to obtain: 


= J |^IF^(l,l')C,f^C|r,.| X (28) 

(jEE ~2 

1 - ' een '^ (1 + - I'D - - I'D) 


where p is the correlation coefficient for the fiducial 
model. We can now calculate the derivatives of the de¬ 
lensed B-mode power spectrum (Eq. 28) with respect to 
the ai parameters. Using the chain rule it can be seen 
that 




dai 


- 2 E 

iGh 


Q^BB,lens,N 




f^K.K~2 

'-'I Pi 


(29) 


where we have rewritten the expression in terms of the 
full-sky derivative with discrete multi¬ 

poles /. We evaluate the lensing derivative in this equa¬ 
tion by differentiating the analytical full-sky expression 
of [i^ for which is correct to leading order in 

C””. The N superscript on the derivative indicates that 
when the beam is large and the E-mode is not signal dom¬ 
inated, the derivative must be evaluated after applying a 
filter (jEE i(jEE,n E-mode power. The derivatives 

of the CIB-lensing cross-power with respect to ai and the 
B-mode power with respect to r are much simpler 


dCf - , 
da. ’ 


dc; 


SB,del 


= Cf^-7r = I) (30) 


dr 


where we have noted that = r x = 1) for 

all r. 

We now turn to determining errors on the measure¬ 
ments. For a measurement of the B-mode power spec¬ 
trum, the errors are determined by the covariance ma¬ 
trix: 


Cov 


BB.BB 

la,lb 


2<5/„ 
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Uy{2la + iy 


BB,res 


■E 


K 
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dCf^ 
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r (31) 

BB.lens 


Cov 


Ll 


dCf^ 
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dc^i^ 


(27) 
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where the second and third lines give the dominant 

S t of the non-Gaussian contribution to the covariance 
. When calculating delensed constraints, we have as¬ 
sumed that the off-diagonal parts of the covariance ma¬ 
trix are reduced by the constant delensing improvement 
factor ao = (Cf'®’'®“)i<ioo/(C'fneglecting 
any scale-dependence. We also neglect the weak depen¬ 
dence of and ao on ai given current constraints. 

We assume that the leasing covariance itself is Gaussian 


KK.KK, 

Cov, 


/sky(2G + 1) 




(32) 


and that an analogous expression holds for E, the un- 
lensed E-mode. We again use first order analytical 
derivatives of the B-mode power spectrum in construct¬ 
ing the off-diagonal part of the B-mode covariance ma¬ 
trix (Eq. 31), unlike in (^ . where full numerical deriva¬ 
tives are evaluated. Higher order corrections can cause 
changes to the off-diagonal covariance, but as we find 
that the omission of the non-Gaussian covariance con¬ 
tributions only causes negligible (much less than 10%) 
changes to our final results, our simple approximation of 
the off-diagonal covariance appears sufficiently accurate 
for our estimates. 

We must also make an estimate of the measurement 
uncertainty on the GIB-lensing cross-spectrum over the 
range of scales that contributes most to the signal-to- 
noise. We consider two cases for the errors on the GIB- 
lensing cross-power, corresponding to current or near¬ 
future measurements: first, a constraint similar to the 
currently available Planck GIB-lensing cross-correlation 
data, second, a substantial improvement on this with the 
errors halved. We will make a simple model of the Planck 
measurement errors as follows. We assume 23 bandpow- 
ers, of width Al = 63.5 (half the width of Planck band- 
powers), with the first bandpower centered at I = 72.5, 
i.e. roughly spanning 60 < Z < 1500 (we maintain our 
cutoff below I — 60). We take the Ui parameters to have 
the same effective bins, but will revisit the choice of ef¬ 
fective bandpower width at the end of this section. The 
magnitude of the errors on each of these bandpowers is 
assumed to be a constant A{2ljCj^/y/2) = 0.02uK.sr, 
which is a good approximation to the the errors on the 
Planck cross-correlation (0, e.g., Table 2 and Figure 
5; optimistically, we only use the statistical part of the 
errors). This corresponds to a ~ 50a measurement of 
the cross-power, with fractional errors on the bins rang¬ 
ing from 6 — 18%. The other, more futuristic, scenario we 
consider corresponds to the same assumed measurements 
but with the errors halved. We assume that any signal co- 
variance between the delensed B-mode measurement and 
the measurement of the lensing-GIB cross-correlation is 
negligible, because the Planck lensing measurements are 
still reconstruction-noise-dominated on the scales most 
important for delensing (and in addition, the areas of 
sky used for the two measurements will often differ). 

With our knowledge of the errors on both observables 
(B-mode power and the GIB-lensing cross power) and 


the dependence of the observables on the parameters 
9p = (r, Ui), we construct the Fisher matrix of Eq. 27 
and obtain the marginalized error on r as follows: 


^marginalized/delensed(^) \/ (-^ ^ 


(33) 


We can now compare two cases. First, a case in which 
we do not delens; in this case, there is no reduction of the 
lensing noise in a measurement of the B-mode power, but 
there is also no delensing-induced degeneracy between 
r and at. Second, a case in which we have delensed; 
this reduces the B-mode noise by up to « 2.2, but also 
introduces a degeneracy of r with the set of uncertain 
GIB parameters Oi. To quantify this tradeoff, we define 
the marginalized improvement factor as the ratio of the 
errors on r for these two cases, now including the effects of 
marginalizing over the uncertain correlation parameters 
Qi in the delensed error: 


Q^marginalized ^0 (^)/f^marginalized/delensed(^) (^^) 

We can plot this ratio as a function of /sky, initially 
assuming that due to foregrounds and systematics, the 
largest scales below = 20 (and hence any tensor B- 
modes from reionization) cannot be measured. 

To gain intuition, we begin by neglecting instrumental 
noise in the B-mode measurement, and simply impose a 
sharp cutoff Z^^ beyond which B-modes are not mea¬ 
sured (though we later investigate more realistic scenar¬ 
ios with different beam and noise specifications). The 
results are shown in Fig. 6, for different ranges of Z which 
contribute to the B-mode power spectrum measurement 
< ^mfx = 100, = 300, Z^f^ = 450) as well as 

two different assumptions for the errors on the measure¬ 
ment of the GIB-lensing correlation (Planck, Planck/2). 

We first discuss our calculation when we restrict 
ourselves to < 100. It can be seen that as 

the area of the survey becomes large, amarginaiized = 
o'o(T’)/o-niarginaiized/deiensed(T’) falls rapidly. As the area is 
increased, the marginalized cr^arginaiized/deiensed(?') must 
therefore become much greater than its value for fixed 
GIB parameters. (Note that ao, the error when not de¬ 
lensing, is unaffected by marginalization). We can under¬ 
stand this as follows. With larger survey area, all con¬ 
straints on r become tighter; with stronger constraints, 
the sensitivity of a constraint on r to cross-power un¬ 
certainty increases, because smaller r can more easily 
be mimicked by uncertain cross-powers as in Fig. 5. As 
large, powerful surveys are more sensitive to cross-power 
uncertainty, marginalization over the uncertain cross¬ 
powers leads to a significant increase of errors on r when 
surveys are large. For a survey larger than ~ 30000 
square degrees, this sensitivity is such that it is no longer 
worth delensing if this comes at the cost of introducing a 
degeneracy of r with uncertain GIB cross-powers. (How¬ 
ever, note that our current delensing weighting is not 
optimized to account for this marginalization.) 

The conclusion that it is not worth GIB delensing 
the largest surveys changes either if the cross-correlation 
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FIG. 6. Delensing improvement factor a on a constraint on 
r, as a function of survey area - now marginalized over the 
uncertain CIB-lensing cross-spectrum. If the CIB-lensing cor¬ 
relation were known perfectly, a would be constant with area, 
as both delensed and non-delensed errors would change by the 
same factor as the area is varied. However, marginalization 
changes this for large surveys, as can be seen from the falling 
blue solid lines (corresponding to Planck-like errors on the 
CIB-lensing correlation) and green dashed lines (correspond¬ 
ing to Planck errors halved). The degradation in delensing 
performance seen is due to a delensing-induced degeneracy 
of r with the CIB-lensing cross-spectrum; for large, powerful 
surveys, this degeneracy significantly increases the delensed 
errors on r upon marginalization. To gain insight, we have 
here neglected instrumental noise, and have simply added a 
sharp cutoff to which the B-mode power is measured. For 
each set of lines, the lower line corresponds to = 100, the 
middle line to = 300 and the upper line to = 450. 
The constant value of a when there is no uncertainty in the 
CIB spectra is also shown as a black semicircular point at 
the top left of the figure. As the resolution of the B-mode 
power spectrum measurement is increased from = 100 to 
= 450, it can be seen that the degeneracy of r with the 
cross-correlation parameters is broken as the CIB-lensing cor¬ 
relation is increasingly self-calibrated by the B-mode power 
spectrum at small scales. 

measurements are greatly increased in accuracy, or if we 
use the constraining power of smaller scales in 
Increasing Zmax to 450, we can see that now the im¬ 
provement upon delensing remains close to 2.2, the value 
we obtained when neglecting uncertainty in the CIB- 
lensing cross-power, even for large surveys. An expla¬ 
nation for this observation is that the degeneracy of the 
cross-correlation parameters and r is being broken by 
measurements of the small-scale B-modes; effectively, the 
higher-Z B-mode power spectrum measurements are de¬ 
termining the CIB-lensing correlation and are fixing the 
amplitude of the low-Z lensing B-mode residual, uncer¬ 
tainty in which would otherwise be degenerate with r. 



FIG. 7. Delensing improvement factor a on a constraint on 
r as a function of survey area, marginalized over the CIB- 
lensing cross-correlation, as in Fig. 6. Instead of the sharp 
ZmaL cutoffs used to gain understanding in Fig. 6, we show 
more realistic results for different noise and beam specifi¬ 
cations {highRes corresponds to a 1.4' beam, lowRes a 30' 
beam). Planck-like errors on the CIB-lensing correlation are 
assumed. For reference, the ideal improvement factors when 
there is no uncertainty on the CIB spectra are indicated with 
solid circles at the left of the plot. Our results again confirm 
that marginalizing over the uncertain CIB spectra does not 
greatly degrade constraints for small surveys or for surveys 
which can resolve small B-mode scales (and can hence self- 
calibrate). For many surveys, CIB delensing thus does not 
require the assumption of a model for the CIB spectra. 


We confirm this by noting that the constraints on any , 
marginalized over r and holding the other a parameters 
at their fiducial values, are tightened significantly as Zmax 
is increased beyond /max = 100 (whereas for /max = 100, 
including the B-mode measurements does not tighten the 
constraints on Oi). It is possible that this self-calibration 
could also have been equivalently achieved via noisy lens¬ 
ing reconstruction using only / = 100 — 450 information 
(though a reconstruction may have insufficient signal-to- 
noise per mode to delens, it may still allow a reasonable 
measurement of the cross-correlation given enough sky 
area). 

As expected, the constraints on r are much less affected 
by marginalizing over uncertain CIB-lensing cross-powers 
if the cross-powers are measured twice as accurately by 
a futuristic lensing survey. 

To make our analysis more realistic, we now repeat 
our calculations for the highRes and lowRes reference ex¬ 
periments introduced earlier, now including instrumental 
noise in our calculation. We show results for both exper¬ 
iments, at two different noise levels, in Fig. 7. We again 
see that for small surveys, marginalizing over the uncer¬ 
tain CIB spectra barely degrades constraints from their 
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ideal values. It can also be seen that the degradation 
with sky area is less severe for high-resolution experi¬ 
ments (blue curves) than for low-resolution ones (green 
curves). This is consistent with the Z^^^-dependence dis¬ 
cussed in the context of Fig. 6, and is due to the fact that 
high-resolution experiments can use the high-Z B-modes 
for better self-calibration. For the lowRes experiment, 
the improvement factor falls, though generally does not 
degrade as much with sky area as when a sharp cutoff 
is applied (as in Fig. 6), because even high-Z B-modes 
that are individually noisy contribute when averaging 
over many modes. Our results confirm that, for a range 
of realistic surveys (in particular, small surveys and high 
resolution surveys), a specific CIB spectrum need not be 
assumed; the spectra can instead be marginalized over 
without a significant loss of delensing performance. 

We will briefly discuss how our results change if we now 
assume that the reionization B-modes at < 20 can 
be measured. In this case, we do not find a substantial 
reduction of delensing performance when marginalizing 
over the CIB spectra, even for very large surveys. We 
explain this as follows: when reionization B-modes can 
be measured, most of the constraining power on r arises 
from the largest scales < 20; on these scales, the 
effects of uncertain CIB spectra are much less important 
(as they are subdominant to the large lensing B-mode 
variance). 

Finally, we briefly comment on our use of effective 
bandpowers as a very simple non-parametric descrip¬ 
tion of the CIB cross-power for forecasting and marginal¬ 
ization purposes (see |59l | for a similar analysis in a differ¬ 
ent context). We focus first on choice of their width. It 
should be noted that to obtain the residual B-mode power 
of Eq. 14, the correlation is integrated over the kernel 
plotted in Fig. 1. Our bandpower widths are chosen to 
be somewhat smaller than the scale on which this kernel 
varies strongly. Still narrower bins are thus unnecessary 
to characterize effects on the residual B-mode power, as 
any sub-bandpower features in the CIB cross-correlation 
will be averaged out by the integration. Of course, since 
most physical processes should lead to smooth features in 
the cross-power, the use of wider bins, different general 
parametrizations that encode smoothness (such as a lim¬ 
ited set of Fourier or polynomial coefficients) or Gaussian 
processes may be well-motivated. However, investigating 
the impact of our choice of bandpower width by doubling 
(and halving) their width, we find that our results are ef¬ 
fectively unchanged. We thus believe that marginalizing 
over our current set of CIB parameters should give a con¬ 
servative, though reasonable, estimate of the degradation 
of constraints on r. 


B. Foregrounds and Non-Gaussianity of the CIB 
Maps 

We have argued in section 3.1 that our forecasts are re¬ 
alistic for many surveys, despite the presence of galactic 


dust foregrounds in the CIB maps. We have also ar¬ 
gued that the raw CIB auto-spectra (including residual 
foreground power) can be measured accurately, so that 
additional foreground power should not bias or degrade 
results if r and the true CIB power spectra are simul¬ 
taneously constrained. However, there are two potential 
limitations of our analysis, which will need to be revisited 
in future work. 

First, we have assumed the CIB maps to be Gaussian, 
whereas some non-Gaussianity is in fact expected, both 
due to intrinsic non-Gaussianity in the CIB field and due 
to residual foregrounds or secondary anisotropies. If we 
assume perfect CMB cleaning, any non-Gaussianity in I 
is not expected to lead to biases, as the mean residual B- 
mode power spectrum is obtained from two powers of the 
GIB field I. Though non-Gaussianity in I could poten¬ 
tially affect the size of the error bars, integration of many 
independent modes is expected to reduce the already- 
moderate non-Gaussianity present in the two-halo-term 
dominated Z = 200 — 500 GIB scales most important for 
delensing. However, this will need to be quantified in 
future work. 

Second, we have assumed that the E and B-mode maps 
have been cleaned successfully using multifrequency data, 
so that any foreground residuals can be neglected. If 
the CMB component separation and cleaning is imper¬ 
fect, additional complications can arise. Though CIB 
anisotropies are expected to be negligible in the large- 
scale CMB polarization (as galaxy polarization directions 
add incoherently), any remaining galactic dust intensity 
features in the CIB may be correlated with polarized 
galactic dust residuals in the CMB E and B maps. In 
principle, this could lead to additional Gaussian and non- 
Gaussian contributions to the delensed B-mode power 
spectrum and its error. While dust residuals are expected 
to be small in the cleaned CMB and CIB maps, and the 
galaxy-dominated low-Z modes (which overlap with the 
scales used to constrain r) can be excluded from the de¬ 
lensing CIB map without significant losses, the magni¬ 
tude of any biases and corrections are currently uncer¬ 
tain. 

Therefore, though simple arguments suggest that the 
impact of foreground residuals and non-Gaussianity on 
the delensing procedure is likely to be moderate, a more 
detailed investigation of these effects, with simulations 
and data, will be required in future work. We note that 
further investigation of some of these potential issues is 
also well-motivated for internal delensing methods. 


V. CONCLUSIONS AND OUTLOOK 

We have examined in detail the methods and perfor¬ 
mance of GIB delensing, including combinations of dif¬ 
ferent frequencies as well as combinations of CIB maps 
with lensing reconstruction. Our results show that CIB 
delensing can remove more than 50% of the lensing 
power and can hence lower the noise of a B-mode mea- 
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surement by a factor of ~ 2.2, and significantly more 
when CIB maps are combined with lensing reconstruc¬ 
tion. Marginalizing over the true CIB spectra, we find 
that uncertainty in the CIB spectra does not degrade de- 
lensing performance significantly for small surveys, and 
only moderately reduces delensing performance for large 
surveys with high angular resolution. This confirms that, 
in most cases, CIB delensing can be performed without 
assuming a model for the CIB spectra. However, further 
investigation into systematic errors from foregrounds and 
non-Gaussianity will be required in future work. 

At noise levels of ~ l^K-arcmin (as targeted by CMB 
stage-IV) or below, internal delensing by reconstructing 
the lensing field from the CMB itself is significantly more 
powerful than CIB delensing. However, for many experi¬ 
ments that are currently planned, CIB delensing appears 
competitive with internal delensing. Even experiments 
such as COrE, LiteBIRD or PRISM [s^, [H, are ex¬ 
pected to have polarization noise levels of 2^K-arcmin or 
higher (especially after foreground cleaning). At these 
noise levels, internal delensing does not perform much 
better than CIB delensing (e.g., in the case of iterative 
delensing with 2/iK-arcmin noise, the improvement fac¬ 
tor is a ~ 2.7 from [s^, whereas CIB delensing gives 
a ~ 1.9). Perhaps most significantly, many current and 
upcoming surveys in the next decade do not have the 
high angular resolution required to delens well internally 
(e.g.. Keck Array, BICEPS, LiteBIRD; so that 

external delensing is the only feasible method; for these 
experiments, CIB delensing is a promising option. 

CIB delensing can significantly strengthen constraints 
on tensor B-mode polarization while requiring only min¬ 
imal assumptions about CIB properties. It is thus a ro¬ 
bust and promising method for enhancing the upcoming 
generation of CMB polarization surveys. 


ACKNOWLEDGMENTS 

We thank R. Keisler, U. Seljak, A. van Engelen, 
O. Dore, J. C. Hill, E. Linder, S. Ferraro, J. Errard and 
J. Dunkley for helpful discussions and/or feedback on the 
paper draft. We are particularly grateful to G. Lagache 
for making the CIB halo model from available. BDS 
was supported by a Fellowship from the Miller Institute 
for Basic Research in Science at the University of Califor¬ 
nia, Berkeley. BDS also acknowledges the hospitality of 
the Beecroft Institute of Particle Astrophysics and Cos¬ 
mology at Oxford University, where part of this work was 
completed. 


Appendix A: Combining Tracers for Delensing 

We now derive the optimal combination of multiple 
fields which trace large-scale structure, such as CIB maps 
measured in different frequency channels (without as¬ 
suming redshift distributions for these fields), in order 


to maximize delensing performance. This optimal com¬ 
bination can be described by writing the tracer J as a 
linear combination of different fields 


/(n) = CiA 


(Al) 


where A are the different fields to be combined, Ci are 
the linear combination coefficients and repeated indices 
are summed over. 

To construct a field which is optimal for delensing 
and cross-correlating, we now solve for the coefficients 
Ci which maximize the correlation coefficient: 


. 2 _ 


kI\2 
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[Ci{Ii X k)Y 


c^Cf 


crc" CY'^CiCjih X Ij) 


We maximize by taking the derivative of In with re¬ 
spect to Ca 
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(A3) 


We note that this equation can be rewritten as a matrix 
equation involving the optimal vector of coefficients or 
c: 


0 = 


2 


k 

c • k 



(A4) 


where we have defined the vector k and the matrix Q as 
follows: 


II 

?: 

(A5) 

Qb = 

(A6) 


It can thus be seen that the correlation coefficient is max¬ 
imized when 


k = Q c X n (A7) 

where n is an arbitrary normalization factor, set by 
Eq. 13. For clarity, we will assume n = 1 for now, and 
will verify that our normalization is correct by checking 
at a later stage that the normalization Eq. 13 is satisfied 
by this choice. 

The optimal coefficients are thus given by 

c = Q-i k (A8) 

Rewriting the original expression for the correlation coef¬ 
ficient in our matrix formalism and inserting the solution, 
the maximal correlation coefficient is given by 


(c • k)^ k • [Q ^ k] 

CY^c ■ [Q c] ^ Cp 


(A9) 
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Inspired by this result we rewrite the weights in terms 
of the correlation coefficients of two li fields, Py , and 
the correlation coefficients of the li fields with the true 
lensing convergence, We can rewrite the matrix Q 
as (no sum) 


The correlation coefficient for this optimal combination 
is 

P — 'y ] PiK.{P )ijPjii (^14) 

b 


^lili 


Qrj = iQ~% = iP~%/4F^l 

(AlO) 


and 

h = P.. X (All) 

We obtain the following weight for a field A (indices 
here are not summed over unless indicated): 


The power spectrum of this linear combination is 


= 1 [(p %3P]k\ [(P ^)abPbK] I X Cf 

yiajb 


PiK (P )ijPjK 


X Cf 


(A15) 


3 



SO that 


^ = E 

i 


y^(P )i3P3K. 

3 



and the cross-spectrum with k takes the exact same 
(A12) value. Hence 

= p^Cd (A16) 

We note that because a = Ci^= 1, the linear com¬ 
bination is already correctly normalized (for our choice 
(A13) of n = 1). 
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